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CN) 1 Context. Within the cosmological concordance model, Cold Dark Matter (CDM) subhalos form the building blocks 

which merge hierarchically to more massive galaxies. This concept requires that even nowadays massive galaxies are 
surrounded by numerous subhalos. Since intergalactic gas is accreted by massive galaxies, observable e.g. as high- 
velocity clouds (HVCs) around the Milky Way, with extremely low metallicities, these can be suggested to represent 
^ the baryonic content of primordial Dark Matter (DM) subhalos. Another possibility of their origin is that they stem 

from disrupted satellite galaxies, but in this case, these gas clouds move unaccompanied by a bound DM structure. 
Aims. HVCs are mainly observed at significant distances from the galactic disk, while those closer to the disk have 
on average lower infall speed and larger but not yet solar element abundances. This can be caused by interactions 
with the hot halo gas: decelerated by ram pressure and metal-enriched by the more abundant gas in the galactic halo. 
Since HVCs are observed with long gas tails and with irregular substructures, numerical models of gas clouds passing 
through hot halo gas are performed aiming at exploring their structure and compare them with observations. If HVCs 
\^ ■ are engulfed by DM subhalos, their gas must leave the DM gravitational potential and reflect this in their dynamics. 

On the other hand, the evolution and survival of pure gas models must be tested to distinguish between DM-dominated 
and DM- free clouds and to allow conclusions on their origin. 

Methods. A series of high-resolution hydrodynamical simulations using the adaptive-mesh refinement code FLASH V3.2 
Q ■ are performed for typical HVC masses, distances, and infall velocities. 

$H ' Results. The models demonstrate that purely baryonic HVCs with low masses are disrupted by ram-pressure stripping 

| and Kelvin-Helmholtz instabilities, while more massive ones survive, losing their initially spherical shape and develop 

, significant substructures including cometary elongations in the column density distribution ("head-tail structure"). On 

the contrary, HVCs with DM subhalos survive with more than 90% of their gas mass still bound and spherically shaped, 
approaching the Galactic disk like bullets. In addition, we find that velocity gradients along the cometary head-tail 
structures does not necessarily offer a possibility to distinguish between DM-dominated and purely gaseous HVCs. 
Comparison of models with observations let us conclude that HVCs are not embedded in a DM substructure and do 
not trace the cosmological subhalo population. 

■ Key words. ISM: clouds - Galaxy: halo - Galaxy: evolution - Cosmology: dark matter 

rn; 

O ' 1. Introduction incompatible with the Galactic rotation ([Blitz et al.lll999t) . 

fN| . This obtrudes the extragalactic origin of these gas clouds, 

t-H ■ In his general investigation of a possible Galactic corona, excluding the wid ely discussed possibili t y of ejecta fro m 

•• ■ |Spitzer| (l95g) was the first who predicted the existence of galactic fountains (|Shapiro fc Fieldfl97ol lBregmanl[l980l) . 
. £h high- velocity clouds (HVCs) . A few years later the system- 
k> . atic search for neutral hydrogen clouds succeeded and E. 

H i R aimond observed the 21 cm line of neutral hydrogen at This isgue leaveg Qpen twQ furthcr fundamental p0S sibil- 

03 . high velocities (|Muller || 1963| ). ities: They are . -q relics of cosmological struc ture formation 

Although the possibilities in observations, the thcorct- ina Cold Dark Matter (CDM) cosmology ([Connors et all 

ical understanding of the interstellar medium (ISM), and |2QQ6l ) and therefore distant "b uilding blocks" f rom galaxy 

the facilities for numerical simulations have increased dras- formation in the Local Group ([Blitz et al.lll999t ) or 2) rem- 

tically, yet half a century later, the formation and evolution nan t s fro m galaxy collisions, f rom tidal disruption of dwarf 

of HVCs are still under controversal discussion. The theo- galaxies ([Putman et all 12004) . or /and from close-by tidal 

rics on their origin do not only vary in small details but lead interactions of satellite galax ies ([Gardiner fc Noguchi|[l996l: 

to completely different pictures of the observed clouds, with iBland-Hawthorn et al.lll998l ). The main difference to dis- 

tremendous consequences for galaxy formation and evolu- tinguish between both sources is their Dark Matter (DM) 

tion. They are approaching the Galactic disk with high ve- content: in case 1) clouds should be strongly dominated by 

locities and, depending on their distance, have to pass the DM, but not in 2). Decisive issues could be derived from 

hot coronal gas of the Milky Way halo. detailed observations in connection with complex numerical 

Their radial velocities lie above or around escape ve- simulations, accounting for their distances, physical state, 

locity of the Milky Way and their rotation component is dynamical structure and chemical abundances. 
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Table 1. Distances and Abundances of HVC and IVC 



Cloud 


Distance 


Abundance 
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(rel. to solar) 




Complex A 
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0.02 - 0.4 
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0.8 - 1.8 
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1 


Lower IV Arch 
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« 1 


1 


LL IV Arch 


0.9 - 1.8 


1.0 ±0.5 


1 


Complex K 


< 6.8 


< 2.0 


1 


Pegasus-Pisces Arch 


< 1.1 


0.54 ± 0.04 


1 


Complex gp 


0.8-4.3 


1-2 


1 



Re ferences. ( l )[Wa kkcr (20(3); (2) lThom etaLI (|2008T l; f3) ICollins et~aTl (|2007l b (4) lLockmarJ (|2003l b (5) Scmbac heTaD IpOO 



(6) iBriins et al. I (2002) 



1.1. Observations 

1.1.1. Distances, All-Sky Surveys, General Properties 

Detailed studies measure distances of a few kpc for the ex- 
tende d HVC complexes (|Thom et all 120061: IWakker et all 
120081). e.g. Comple x C with a distance of d = 10 ± 2.5 kpc 
( Thorn et all 120081). For t he massive Complex H a dis- 
tance of 27 kpc (lLockmanll2003l ) and a mass of 2 x 10 7 M Q 
(|Simon et al.ll2006| ) is derived. Beside such large complexes 
which cover around 40 % of the sky for a re solution limit of 
N H i > 7 x 10 17 cm" 2 (jMurphv et al.lll995l) . there have also 
been found isol ated, compact gas clou dlets with a small 
angular extent (jBraun fc Burtonl Il999f l. The first all-sky 
study of this sub-population, called Compact High Velocity 



de Heii et al. (2002a 



de Heii et all (|2002bD 



Clouds (CHVCs), was presented by 
as a combination of data extracted by 
from the Leiden Dwingeloo Survey for the northern hemi 
sphere and fr om the HI Parkes Al l Sky Survey (HIPASS), 
extracted by IPutman et al.l (|2002h for the southern hemi- 
sphere. Due to the compactness of CHVCs, the lack of dis- 
tance indicators, and the poor knowledge about their physi- 
cal environment, direct distance measurements are difficult. 

Observations of M31 have shown recently that HVC- 
analogs arc clustered around th eir host galaxy with a max- 
imu m distance between 50 kp c (jThilker et al.ll2004|) and 60 
kpc ( Westmci cr et al.ll2005bh . Most recently. iRichter et all 
(|201ll ) have studied Ca II absorbers at low redshift and pro- 
pose that the characteristic radial extent of HVC analogs 
with log N (HI) > 17.4 in galaxies with z < 0.5 is Rhvc ~ 
55 kpc. 

Another remarkable property of HVCs is a 2-phasc 
structure with a cold, dense core and a warmer, 
more tenuous enve l ope, w hich was already discussed by 
IWakker fc Schwarzl (|l99lh for clumps in the larger com- 
plexes and is still supported by more recent observa - 
tions of smaller cloudlets as in iBen Bekhti et al.l (|2006l) . 
ICiovanelli et al.l ()2010D presented ALFALFA (Arecibo 
Legacy Fast ALFA) observations of ultra-compact high 
velocity clouds (UCHVCs) at assumed extragalactic dis- 
tances. The derived masses of the measures UCHVCs at a 
distance of 1 Mpc are below 1O 6 M , which explains the 
non-d etection by previous HI surveys, such as lPisano et al.l 
(|2007h . 



1.1.2. Dark Matter Minihalo and Ram Pressure Stripping: 

If the CHVCs permeate the intergalactic space their ori- 
gin allows two possibilities: they are disrupted remnants 
from galaxy collisions and formed in the tidal tails or they 
represent the relics of structure formation in a CDM cos- 
mology. The latter possibility would provide the requested 
link to the postulated DM building blocks of hierarchical 
galaxy formation which means that this baryonic cloud 
mass should be harbored by a so-called DM mini-halo. 
iNichols fc Bland-Hawt horn ( 200^) have investigated the or- 
bit of the Smith Cloud under the consideration that the 
cloud is confined by a DM halo. 

Detailed observations of HVCs are showing a head-tail 
structure in the hydrogen column density di stribution in 
both over-densities within the large complexes (|Bruns et al.l 
20001) and s i ngle CHV Cs, like presented for e xample by 



Briins et all (|200lD and lBen Bekhti et al.l (|2006l) . In agree- 



ment with theoretical and numerical studies this structural 
properties are unambiguously indicative for gas clouds mov- 
ing through a surrounding plasma. It is expected that if 
the force due to ram pressure exceeds the gravitational 
restoring force that keeps the cloud material via self-gravity 
toget her, material is stripped off the cloud (jGunn fc Gottl 
I1972D . This leads to a head-tail structure of the cloud, which 
is observed with radio observations of the 21cm line of 
neutral hydrogen (IBriins et al.l 120001: IPutman et "all 120021; 
IWestmeier et al.l l2005al) . Furthermore, the mixing of hot 
halo and cool cloudy gas might be resp onsible for th e ob- 
served ionization stages in the hot halo (|Spitzerlll99(l . 

Observationally derived structure parameters allow the 
conclusion that CHVCs are bound by a DM halo of signifi- 
cant mass (|Braun fc Bu rton 200(|) . While the gravitational 
potential in the vicinity of the HVC is dominated by the 
combined gravitational potentials of the cloud itself and 
the host galaxy, a DM halo around the CHVC enlarges the 
total mass of the cloud-halo system and therefore also the 
volume in which the cloud potential exceeds the potential 
of the major galaxy, so that HVC material is gravitationally 
stronger bound to the cloud. 

In addition, whether HVCs represent the baryonic con- 
tent of cosmological DM sub-halos or not, should become 
visible during the ram pressure stripping (RPS) action on 
their passage through the hot halo gas. The gravitational 
acceleration in the vicinity of the HVC is dominated by 
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Table 2. Initial parameters of the high- velocity clouds in different simulations. First column: initial setup no. for every 
run; col. 2: relative HVC velocity v w ; col. 3-5: HVC temperature T gasi radius R gas in pc, and mass M gas in Mq) DM 
subhalo scale radius of tq (col. 6) in pc and mass Mum (col. 7) in M Q ; col. 8: baryon-to-total mass ratio M gas / M to u 
maximum refinement level of the adaptive grid (Max.Rcf.) in col. 9 and effective resolution (Res.) in col. 10. In the last 
column (Run Nr.) a unique number is assigned to every simulation run. 



Setup 


V w 


Tgas 


Rgas 


M gaa 


ro 


Mdm 


Mgas/Mtot 


Max.Ref. 
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r wi 
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5 
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6 
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- 


7 


5 


VI 
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127 


4.40 x 10 5 
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0.08 


7 


5 
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4.41) X 1U 
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0.03 


6 


25 
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4.40 x 10 
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1.63 x 10 7 


0.03 


7 


5 


IX 
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4.40 x 10 s 


- 


~ 


- 


6 


18 


X 
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1000 
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4.40 x 10 5 


250 


1.63 x 10 7 


0.03 


6 
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XI 


3 
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376 


3.10 x 10° 








4 


101 


XII 




200 


2000 
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3.10 x 10 6 








5 


50 


XIII 




200 


2000 


376 


3.10 x 10 6 








6 


25 
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376 


3.10 x 10 6 
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3.10 x 10 6 
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4.93 x 10 6 


0.39 
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50 
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25 
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50 
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the self-gravity of the cloud and the external gravitational 
potential of its host galaxy. At first, the cloud gas is con- 
fined within the equipotential of cloud and galaxy poten- 
tial. Once the stripped-off material runs over this region, it 
gains an additional acceleration in comparison to the ma- 
terial that is stripped off but still gravitationally bound 
to the cloud. This effect is expected to be hardly notice- 
able if the HVCs total mass is only determined by its gas 
mass, because the potential of the massive host galaxy mass 
would widely dominate the self-gravity of the cloud. If the 
cloud is embedded in a DM subhalo the total mass of the 
HVC and therefore its self-gravity is much higher and could 
make this effect visible in velocity-position space as a sharp 
change in the velocity gradient along the head-tail structure 
of the cloud, as observation s obtrude (jBruns et al.l l2000l 
[20011 iBen Bekhti et aLll2006h . 

1.2. Numerical Simulations 
1.2.1. Ram Pressure Stripping: 

iQuilis fe Moord (|200lh have studied by means of hydro- 
dynamical simulations the behavior of CHVCs passing- 
through the dilute hot halo gas of a typical disk galaxy at 
large distances. The main purpose of their numerical mod- 
els was to compare the emerging head-tail structure with 
observations, in order to draw conclusions on the necessity 
of a DM component and on a possible lower density limit 
for the surrounding intergalactic medium (IGM). 

They found head-tail structures in environments with a 
density higher than 10 -4 cm -3 , whose lifetimes in simula- 



tions without a DM halo are very short, w lOMyrs only, 
while the tails in their DM dominated clouds survive for 
« 1 Gyr. 

Although the authors have simulated a grid of HVC 
models, it is worth to point out, that there is no compari- 
son of the same cloud with and without DM, like we discuss 
in this work, but the authors study two different scenarios: 
HVCs as close-by structures at a distance of 1-10 kpc and 
masses between 10 and 100 Mq vs. very distant DM dom- 
inated clouds at a distance of 300 kpc, a radius of 1 kpc 
and a mass of more than 10 7 M Q , co ntrary to the typica l 
maximum distance of around 55 kpc ([Richter et alj|20l"ll) . 

A snapshot of the DM dominated cloud modeled over 
300 Myrs (Fig. 2 of IQuilis fc Moord (|200lll ) with a con- 
stant velocity of 200 km/s travels over a distance of 60 
kpc through the galactic halo, i.e. must have faced changes 
in the ambient densities, a tidal field and an acceleration 
by the host galaxy. Due to the lack of self-gravity, low- 
mass DM-free gas clouds are disrupted ver y quickly so that 
the e volution is only shown after 3 Myrs (jQuilis fc Moord 
|2001L Fig. 3). The DM-free clouds in the simulation runs 
presented in this work can easily survive for 100 Myrs, not 
only because of the included self-gravity but also due to 
the much higher masses and distances, in agreement with 
recent observations. 

iHeitsch fc Putmanl (|2009T) have shown with grid-based 
hydrodynamic simulations that smaller clouds with HI 
masses < 10 4 5 Mq will lose their HI content during their 
paths of 10 kpc at maximum (for typical relative velocities 
and halo densities). They propose that the stripped ma- 
terial which is usually gas in the warm and ionized phase 
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can contribute to the exte nded layer of warm, ionized gas 
at lo wer galactic heights (jRevnoldsl Il993t iGaensler et al.l 
120081) and later on rccool and form smaller HI cloudlets 
which are classified as low or intermediate velocity clouds 
(LVC, IVCs). 



Table 3. The initial distribution of the cloud follows in all 
setups a hydrostatic equilibrium, which is derived by Eq. [2 
where M(r) includes either only the gas mass or in addition 
to that also the enclosed DM mass following Eq. [3] and S] 
with a scale radius of ro = 150 pc or tq = 250 pc. 



1.2.2. Heat Conduction: 

Due to the large temperature gradient between the cloud 
and the surrounding halo gas, saturated heat conduction 
should be an important process of energy transport be- 
tween the gas phases, leading to the classical two-phase 
structure. Furthermore, heat conduction stabilizes clouds 
with head-tail structure and is extendi ng their lifetime 
(IVieser fe Henslerl l2007bl) . In addition, |V ieser fe Hensler! 



2007al) have shown that in models which combine heating, 
cooling and saturated heat conduction with self-gravity, 
condensation of halo gas on the cloud surface can domi- 
nate evaporation and lead to a net accretion of ambient 
gas. 

Fo r the warm cosmic rain theory by iHeitsch fe Putmanl 
(l2009h this means that the stripped off material is not only 
warmer and higher ionized than the cold core of HVCs, 
but would have other significant consequences: On aver- 
age HVCs contain much smaller element abundances than 
the Milky Way ISM (see Tab. [1]) what clearly advocates 
for their primordial origin. HVCs affected by heat conduc- 
tion and stripping in the presumably metal-rich halo gas 
are expected to have also a higher metallicity because the 
stripped cloud surface consists not only of cloud gas but also 
of accreted halo gas. If part of the LVC and IVCs are the 
decelerated survivors of disrupted or ablated HVCs they 
should have a higher metallicity than the less processed 
HVCs, which is observationally confirmed (Tab. [1]). Since 
IVCs/LVCs rain down with velocities below escape velocity, 
an origin from the galactic fountain cannot be ruled out. 
The cores of HVCs, which survive ram-pressure stripping, 
can be decelerated by either the drag of the hot halo gas 
or other processes like the buffering of plane-parallel galac- 
tic mag netic field lines whic h have to couple to the HVC 
plasma (|Zimmer et al.|[l997f l. This effect would contrib ute 
to the heating of the halo gas (jjelinek fc Henslerll201ll ). 

In a rece nt study about the origin of the HVCs 

iBinnev et al.l (|2009[ ) suggested that they can be formed by 
thermal instability in the warm-hot Galactic halo medium. 
For Milky Way sized halos, however, they found that this 
happens only in regions outside of 100 kpc and for an al- 
most perfectly fiat entropy profile, when thermal instabil- 
ities are not damped by a combination of buoyancy and 
thermal conduction. They concluded that it is therefore a 
rather unlikely formation scenario for HVCs. 

1.3. General Questions 

Since differences in the ram pressures stripped gas and in 
the head-tail structure are not yet explored with respect to 
their possible origins, here we present numerical models of 
HVCs with and without a stabilizing DM subhalo based on 
the following questions: Is it possible for the cloud with DM 
to keep the gas longer bound as without? Does a DM halo 
suppress hydrodynamical instabilities? Is a cloud without 
a DM halo already disrupted and does it lose most of its 
gas due to the ram pressure produced by its supersonic mo- 
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tion (up to more than 400 km s _1 ) through the hot gaseous 
halo of its host galaxy or in the outskirts of the dilute halo 
gas? Can we distinguish from the observed velocity - po- 
sition distribution of the head-tail structure whether the 
observed cloud contains DM? If they are embedded in a 
DM halo, are HVCs the low mass extensions of gas rich 
subhalos with (ultra faint) dwarf galaxies as their higher 
mass analogs? Moreover, if they are DM-free and a signif- 
icant fractio n of the few observed satellite galaxies have a 
tidal origin ([Bournaud fc Dud [2006; O kazaki fc Taniguchil 
I2000D : where are all the DM satellites we see in large cos- 
mological simulations? 
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Fig. 1. Initial density (top panel) and pressure (lower 
panel) distribution of the model setups described in Tab.O 



4 



S. Plockinger and G. Hensler: Do HVCs trace the DM subhalo population? 



1.4. Paper Structure 

This paper has the following structuring: First the numer- 
ical method and the model setup of both the host galaxy 
and the HVCs are presented. In Sect. 12.41 we discuss the 
size of the Roche-lobes in this 2-body system and the con- 
sequences for the gas that is lost by RPS. 

The effect of RPS that leads to the observed cometary 
appearance of the clouds is discussed in Sect. 12. 51 while 
further possible mass-loss procedures are mentioned in 
Sect. O 

Because HVCs do not necessarily fall directly towards 
the observer nor radially into the center of the Milky 
Way we take a close look on the line-of-sight problem in 
Sect. 13.21 Not only the column density distribution varies 
when changing the inclination between object and observer 
but the velocity measurements can cover only the radial ve- 
locity. Therefore the whole phase-space distribution of the 
observed clouds is highly sensitive to line-of-sight effects. 
For all performed simulations we processed the data in a 
way that we simulate different line of sight directions to 
look for signatures that are dependent on the line of sight. 

Finally, the dependency of the results on parameters 
like the initial velocity (see Sect. 13.3. ID and the DM mass 
(see Sect. I3.3."2|) and shape (Sect. I3~3U|) are discussed. 

We have included not only a self-gravitating gas cloud 
as well as a static gravitational potential that represents a 
DM mini-halo but also an external gravitational potential 
which gives us the opportunity to study the equipotential 
lines of the combined potential of the galaxy and the HVC. 

The initial distance between the HVC and the external 
mass was chosen to be 50 kpc, which leads to cloud masses 
in the order of 10 5 to 10 6 solar masses. Outside of 60 kpc 
the HVCs are either fully ionized and therefore not visible 
in neutral hydrogen or HVCs origin from an interaction sce- 
nario like the Magellanic Stream which would also explain 
the typical maximum distance of around 50 - 60 kpc. 



2. Models 

2.1. Numerical Methods 

In order to resolve hydrodynamical instabilities and tur- 
bulence in the stripped-off gas, but not to waste compu- 
tational time and spatial resolution by the dynamics out- 
side of the clouds, an adaptive-mesh algorithm is a proper 
choice. 

For all simulations in thi s project the hydro dynamic 
code FLASH3.2 developed bv lFrvxell et~aTI (120001 ) is used, 
which alrea dy prov i des a fast multigrid Poisson solver which 
is based on lRickerl (|2008f ) and adapted for the FLASH grid 
structure by the Flash Center. 



2.2. Initial Setup 

The initial setup for every run is an isothermal cloud in hy- 
drostatic equilibrium with te mperature T nns . radius R gas 
and mass M gas embedded in a IBurkertl (|1995f ) halo with a 
scale radius of r and a mass Mum- The baryon-to-total 
mass ratio for the cloud is given by M gas /M to t and the rel- 
ative velocity between the cloud and the ambient medium 
is given as v w . 

The simulation of an object moving towards the cen- 
ter of an external gravitational potential representing the 



host galaxy, is performed in a way, so that the rest frame 
remains at the cloud's mass center. This leads to a wind- 
tunnel setup with an increasing wind velocity, because the 
gravitational attraction of the host galaxy increases as the 
cloud approaches. Since observations do not show any star 
formation in HVCs, neither nowadays nor at any point in 
the past, the simulated cloud should not collapse due to 
self-gravity during the simulation time, but also not expand 
and dissolve at the simulated heights because we focus on 
the compact, isolated sub-class of the HVCs. Therefore we 
assume an initial hydrostatic equilibrium: 



dP 

dr 



-p{r) 



<9$(r) 
dr 



(1) 



where <I>(r) is the gravitational potential of the cloud and 
since the cloud is spherically symmetric in the beginning 
the pressure distribution follows: 



dP _ d GMjr) 

dr dr r 



(2) 



where p(r) is the gas density and M{r) is the enclosed mass 
within the radius r. 

We intend to compare clouds that arc embedded in a 
DM halo to clouds that only consist of baryonic matter. In 
these two cases the initial pressure and therefore density 
distribution are different since the enclosed mass contains 
both dark and baryonic matter. Subsequently other param- 
eters as the cloud radius and the cloud mass differ between 
the cases with and without a DM subhalo. 

To avoid effects that arc due to the different initial dis- 
tributions we set up clouds that are in hydrostatic equilib- 
rium with the gaseous part only (Setup 2 and Setup 3) as 
well as clouds that are stable under both the baryonic and 
the DM mass (Setup 2-D, Setup 2-DD and Setup 3-D) 
and performed simulations with and without a DM subhalo 
for each setup (see Tab. [2] and Fig. [T]for details). 

Stability studies without wind are performed to ensure 
that the timescale of the changes due to the deviation from 
the hydrostatic equilibrium is negligible, i.e. much larger 
with respect to the ram pressure timescale. 

The cloud is embedded in a homogeneous medium with 
a temperature of Thim — 10 6 K and a density of Phim = 
4 x 10~ 4 H enr 3 which represents the hot galactic halo. 
While observationally determ ined values for the oute r halo 
of the Milky Way are rare, iQuilis k Moord (|200lh have 
shown that head-tail features only appear if the external 
medium has a density larger than 10 -4 cm~ 3 and a relative 
velocity of 200 km s _1 . 

The cloud's pressure is integrated following Eq. [5] down 
to a radius where the pressure of the cloud equals the pres- 
sure of the ambient medium. 

Since our study is aiming at exploiting whether HVCs 
would represent the extension of DM-dominatcd baryon as- 
semblies from dwarf galaxies to the lower masses it seems 
reasonable to choose also as DM den s ity di stribution the 
empirically derived profile by IBurkertl (jl995h : 



Pdm = 



(r + r )(r 2 + r 2 ) 
with the empirical correlation between po and of 

-2/3 



(3) 



po = 4.5 x I0- 2 L^-j M oP c- 3 



(4) 
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Fig. 2. Column density distribution of the gas cloud overlaid with equipotential lines for Setup 2-DD (sec Tab. [5] for its 
parameters) at a maximum refinement level of 7 (effective spatial resolution: 5 pc): without DM halo (left panels, Run 
VI) and with DM halo of scale radius ro = 250 pc (right panels, Run IX). The contour lines of the greyscale Hi column 
density distribution range between \og(Nni) = 18.8 ... 21.2 cm~ 2 in steps of 0.3 dex. Overlaid are the equipotential lines 
of the total gravitational potential by the gaseous cloud, the DM subhalo of the cloud (if present), and the external galaxy 
(see text) in steps of 2.5 x 10 11 cm 2 g~ 2 . In red, the equipotential line of the innermost Lagrangian point is depicted. Top 
panels show the initial distributions, bottom panels after a simulation time of 100 Myrs. (Step-like contour lines reflect 
the coarser grid further away from the cloud.) 



In contrast to th e mostly referred DM profile, the so- 
called NFW profile (|Navarro et al.lll997t) . which was ini- 
tially derived from large-scale cosmological simulations, 
with the characteristic cu spy center (see also iMoore et al.l 
119991 : iNavarro et aLll2004D . the observed rotation curves of 
dwarf galaxies show nearly solid body rotation, which is 
prefer ably fitted with a central core in the DM distribution 
(e. g !Burkert|ll995l: Ide Blok fc Bosmal 120021: iGentile et ai] 
120051: Ide Blokll2005l) . The resulting DM halos in more recent 
high resolution cosmological simulations (|Navarro et al.l 
l2010f) support central cores that can also be derived from a 
global D M-radius relation f rom spirals to dwarf spheroidal 
galaxies (IWalker et al.ll2010D . For further discussions on this 
topic see Sect. 13.3.31 

We decided to run all our simulations up to a given 
simulation time of 100 Myrs, so that we can neglect the 
change in the ambient pressure, while the effects of the 
time-dependent external potential arc well described in our 
simulations. The distance of the CHVC on radial trajectory 
is in our case between 30 kpc and 50 kpc. We expect that for 
clouds closer than that, the simplification of a radial infall 
is not valid anymore and an eccentric orbit close to the disk 
would lead to tidal stretching of the clouds as observed in 
the large complexes closer to the Galactic disk. 

In order to cover a large range of parameter space but 
in a reasonable expense of computational time and to per- 
ceive clear differences between DM-dominated and pure gas 
clouds we use only two different scale radii (ro = 150 pc and 
ro = 250 pc) for the DM profile. In combination with a sam- 
ple of baryonic masses we obtain setups with very different 
baryon-to-total mass ratios = Alt/Mtot both including 
very low as well as very high baryon contents. 



The larger DM halo with r = 250 pc has a total DM 
mass which is more than 3 times higher than that of the 
smaller DM halo with ro = 150 pc. Due to the different 
scale radii and the core shape of the potential, the ratio 
of the enclosed mass of the central part does not represent 
the full mass ratio of the total halos. But within 400 pc, a 
typical extent of the head-tail structure, the larger halo has 
already 70% more DM mass enclosed than the smaller one, 
leading also to a difference in the gradient of the DM po- 
tential of 70%. If the simulations were very sensitive to the 
scale radius of the surrounding DM halo, we would expect 
a significant difference in the results. 

iHoeft et al.l (|2006f ) studied structure formation by 
means of hydrodynamic simulations for different spatial res- 
olutions and find that low-mass haloes below a critical to- 
tal mass of M c (z = 0) = 6.5 x 10 9 /i _1 M & can be baryon 
poor with fb as low as 0.01. We cover the full range from 
baryon poor haloes over haloes with the universal baryon- 
to-total mass ratio d erived by WMAP3 data (ft, w 0.17, 
iKomatsu et al.|[2011[ ) to systems with ft, as high as 0.39 
(Setup 3) to investigate the influence of a DM subhalo in a 
RPS scenario. The latter anticipates that the outermost 
part of the DM subhalo is already distorted by the Galactic 
tidal force. 

2.3. Parameter variants 

For a better comparison between the simulation runs with 
and without DM content for every mass range runs with 
different stability setups are performed. One setup repre- 
sents an initial distribution in hydrostatic equilibrium for 
a self-gravitating gas cloud (setup 2 and 3), while another 
setup is stable in a static run that also includes DM with 
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Fig. 3. HVC model with an initial gas mass of M gas = 4.4 x 10 5 M Q (Setup 2-DD): color-coded pressure distribution 
(logarithmic scale [dyne cm -2 ]) through the midplane (z=0) of the 3D simulation box overlaid with Hi column density 
contours between log(Njji) = 18.5 ... 21.5 cm -2 in steps of 0.5 dex, 100 Myrs after the onset. Upper left panel: cloud 
without DM (Run VI), upper right panel: cloud embedded in an NFW halo (see Sec. I3.3.3j) . lower panels: clouds with 
Burkert halos with ro = 150 pc (Run VII, left panel) and ro = 250 pc (Run IX, right panel), respectively. The red contour 
depicts the (Eki n + E{) = —E pot boundary and encloses all material, both from the cloud and from the hot halo gas, that 
is gravitationally bound to the HVC. The velocity vectors indicate a maximum velocity of 215 kms -1 and their lengths 
are linearly scaled. 



a mass of 4.93 x 10 6 M Q and a scale radius of 150 pc for 
setups 2-D and 3-D and with a DM mass of 1.63 x 1O 7 M 
and a scale radius of 250 pc for setup 2-DD (see Tab. [3]). 
For all these initial distributions we performed runs with 
and without a DM halo. This leads to configurations that 
are not in perfect hydrostatic equilibrium (see Sec. 12.11 for 
details). In Table [2] a summary of the used parameters can 
be found. 

2.4. External gravitational potential 

Obs ervations of indivi d ual C HVCs by iBruns et all (|2001[ ) 
and iBen Bekhti et all (|2006[ ) show not only a head-tail 
structure in the column density distribution but also a clear 
change in the velocity gradient along the elongated density 
contours. 

Since the gravitational potential in the vicinity of the 
HVCs is dominated by its own mass and the external grav- 
itational held of its host galaxy, this system can be consid- 
ered as a two body interaction. Therefore the area within 
the equipotential line through the Lagrangian point 1 (in 
binary systems called "Roche lobe" ) determines the region 
in which the cloud material is gravitationally bound to the 
HVC. Ram pressure leads to an elongation of the gas den- 
sity profile in wind direction and as soon as gas has crossed 
the Roche lobe it becomes unbound from the HVC and 



should therefore be exposed to an additional acceleration 
visible in the immediate change of the velocity gradient. 
Since the size of the Roche lobe is dependent on the mass 
of the HVC, the observed position of the velocity gradi- 
ent change can help to distinguish whether HVCs contain 
a significant fraction of DM or not. 

In our simulations the cloud will start at an initial dis- 
tance of 50 kpc to an external point mass of Mq a i = 
10 11 Mq which represents the host galaxy. IXue et all (2008) 
have derived the Milky Way DM halo mass from the kine- 
matics of blue horizontal branch stars in the galactic halo. 
They fitted both an unaltered NFW profile as well as an 
adiabatically contracted NFW halo and find virial masses 
between M vir = [0.82±°; 2 ^ ... 1.2lt°;*°] x 10 12 M Q with 



virial radii in the range of r V i 



[258 



+20 
-21 



29311*] kpc 



We compared the gravitational acceleration in the cloud 
region of the external point mass with M Ga i = 10 11 M to 
the NFW profile represented by [M mr = 0.82 x 10 12 M Q , 
r vir = 258 kpc, c = 12.2] and found that the relative dif- 
ference is as small as 1-3 % whi ch is much less t han the 
difference between the models bv lXue et all (|2008l) . 

The equipotential lines together with the column den- 
sity distribution of the isothermal gas distribution initially 
as well as after 100 Myrs of simulation time arc shown in 
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Fig. 4. As Fig. [3] The initial gas mass is M gas = 3.1 x 
10 6 M Q (Setup 3). Upper panel: DM-frec cloud (Run XIV); 
lower panel: cloud with DM subhalo of tq = 150pc (Run 
XVII). 



Fig. [U for Setup 2-DD with (Run VIII and IX) and without 
DM (Run V). 

In all setups the gaseous parts of the clouds exceed the 
smallest Lagrangian equipotential surface in at least one di- 
rection. In the setups without the additional mass of a DM 
halo, the LI surface lies completely within the cloud. Note 
that the Roche lobes indicated in Fig. [2] become smaller 
during time not because of the mass loss but because the 
cloud approaches the galaxy mass with an initial velocity 
of 200 km s _1 and reaches a distance of around 30 kpc after 
the simulation time of 100 Myrs. 

Including an external mass in our simulations does not 
only effect the region where the cloud material is bound 
to the objects, but also influences the motion of the HVC. 
Due to the gravitation by a host galax y the cloud is accel- 
erated towards the external point mass. Beni amm fc Danlvl 
(|l997f) have shown that clouds reach a terminal velocity on 
their way to the galactic disk, where the acceleration due to 
gravitational force is equal to the drag force caused by the 
different densities of the hydrostatic layers of the galactic 
disk. Their model is valid for distances up to 10 kpc and 
they did not take the change in the area which is exposed to 
drag into account. In our setup the cloud travels in galac- 
tic heights between 30 kpc (after 100 Myrs of simulation 
time) and 50 kpc (initial distance), where the change of 
the hot halo gas density is negligible. The stronger gravita- 
tional acceleration at smaller distances and in many cases 
the shrinkage of the cloud surface that is exposed to drag 
due to RPS, together with a constant ambient medium does 



no t lead to a terminal veloc ity in our models as assumed 
bv lBeniamin fc Danlvl (|1997[ ). 



2.5. Ram-pressure stripping 

We calculated the analytical stripping radius R s trip as the 
radius where the force due to ram pressure F rarn equals the 
gravitational restoring force F grav . For radii larger than 
Rstrip the gravitational energy of the cloud mass is not 
large enough to keep the gas bound so that it is stripped 
off due to th e relat ive motion of the hot ambient gas. 
iGunn fc Gottl (|1972l ) derived R s trip for the instantaneous 
RPS of a gaseous disk moving face-on through the intra- 
cluster me dium (ICM). IMcCarthv et all (|2008l) modified 
the orig inal IGunn fc Gottl (119721) formulation to determine 
the stripping radius for spherically symmetric density dis- 
tributions. Their stripping condition then reads analogously 



max 

(R) x £ sos (i?) 



where p am b here is the density of the surrounding halo 
material, v re i is the relative velocity between the cloud and 
the ambient medium, g ma x{R) is the maximum restoring 
acceleration in wind direction and T, gas (R) is the surface 
density of the gas projected in wind direction. R is the dis- 
tance from the cloud center perpendicular to the direction 
of the ram pressure. Fo r a schematic diagram see Fig. 3 of 



IMcCarthv et al.l ((200 



To calculate R s tnp, where 



Pamb X ^rel Qmaxi^Rstrip) X ^gasi^^strip) 



(6) 



at any time, we use the 3D data cube and drill bars out of 
the cloud in the streaming direction of the halo gas with a 
cross section of one grid cell at the highest resolution and 
calculate the right-hand side directly from the gravitational 
potential and the gas density that is already stored in each 
grid cell. The ram pressure (left-hand side) would be con- 
stant for a standard windtunnel setup. In our case p am b 
stays constant over the simulation time but the relative ve- 
locity is increasing because of the gravitational attraction 
of the external mass. In our setups the velocity increases 
by about 6% in 100 Myrs leading to an increase in the 
ram pressure by « 12.4%. Therefore the whole condition 
for RPS shown in Eq. [6] is not constant in time because 
both the velocity v re i, effecting the ram pressure term as 
well as the density distribution of the cloud changes, lead- 
ing to a different column density distribution T, gas (R) and 
subsequently to a different gravitational restoring force. 

2.6. Numerical Resolution 

In numerical simulations reliable results should not depend 
on the spatial resolution but only on the included physics. 
To make sure that the results are valid for our setup we 
checked for convergence in performing the same simulation 
run for different levels of refinement. We expect differences 
as well in the apparent substructure of the cloud as in the 
overall dynamical properties for models with low refine- 
ment levels, compared to runs with sufficiently small grid 
sizes. Especially, mass loss procedures as RPS and Kelvin- 
Helmholtz instabilities are highly sensitive to the physical 
size of a single cell. In addition the equipotential volume 
within the cloud, i.e. where the cloud gas is gravitationally 
bound, is getting smaller when the cloud is approaching 
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Fig. 5. Projected velocity-space distribution of Setup 3 without a DM halo at a maximum refinement level of 7. Top left 
panel: contour lines of the column density distribution (in z direction) at log(-/V# /) = 18.8 ... 21.2 Hcm~ 2 steps of 0.6 dex 
100 Myrs after the simulation onset. Both the x- and the y-coordinates correspond to the cartesian grid of the 3D data 
cube. Radial velocities measurable along the line-of-sight (dashed line) in different directions of angle a = deg, 15 deg, 
and 30 deg with respect to the x-direction are shown in the 3 frames from bottom left (a = deg) , through top right 
(a = 15 deg) to bottom right (a = 30 deg). Along the abscissa, projected distances perpendicular to the line-of-sight 
(dotted line) through the mass center (x=0) are shown. All cells with densities n > 0.01 H cm~ 3 are included and re- 
arranged to represent the column density distribution in the observed position-velocity space (leading to the bubbled 
contour structure). The contour lines of the line-of-sight column densities indicate log(Njji) = 18, 19, 20 and 21 H cm -2 . 
The measured maximum velocities differ accordingly from 215 (for a = Odeg) to 188 kms -1 (a = 30 deg). 



the external center-of-mass and is therefore likely to be not 
properly resolved in low resolution runs. 

In Fig. [5] the evolution of both the cloud radius and 
the stripping radius is shown for different levels of refine- 
ment. For the low resolution runs with 64 3 cells (Ref 4) 
and 128 3 cells (Ref 5) the calculated stripping radius suf- 
fers from the poorly resolved gravitational potential within 
the cloud. Therefore the maximum restoring acceleration 
9max is underestimated. For higher resolution runs with 
256 3 cells (Ref 6) and 512 3 cells (Ref 7) results are con- 
verging. Comparing the analytically determined stripping 
radius to the cloud radius for the same resolution, it is 
evident that the evolution of the cloud radius follows the 
changes in the stripping radius. Therefore, dynamical strip- 
ping due to ram pressure is the dominant mass loss proce- 
dure in the considered scenario. 



2. 7. Further mass-loss processes 

Even though we performed highly refined simulations with 
a resolution as low as 5 pc, it is not easily possible to take 
processes into account that are happening on sub-grid scale. 
In this section, we discuss the relative importance of other 
than the above described mass-loss processes, as namely 
Kelvin-Helmholtz (KH) instability (Sec. l2T7TTj) . Rayleigh- 
Taylor (RT) instability (Sec. 12.7.20 . mass loss due to in- 
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details see text). 



trinsic velocity dispersion (Sec. I2.7.3[) and tidal disruption 
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2.7.1. Kelvin Helmholtz instability 

KH i nstability (for detailed analysis, see IChandrasekharl 
Il96l!) can occur when a velocity shear is present within 
a continuous fluid or, when there is sufficient velocity dif- 
ference across the interface between two fluids. It can be 
an important process that leads to a high mass loss due to 
so-called KH stripping of cloud material. Stability analysis 
of the KH instability shows a condition for the growth of 
perturbation in dependence on the wavenumber k = 2tt/X. 
Perturbations with a wavenumber of k < k m i n are damped 
while perturbations with k > k m i n are enhanced (see 
Eq. [7]). The critical value k rn i n can be derived from 



gz(ai - ag) 



(7) 



with a, = Pi/ (pi + p 2 ) (i=l,2), where g z is the gravi- 
tational acceleration perpendicular to the relative motion, 
pi and P2 are the densities e.g. of cloud and hot gas at 
the edge of the cloud, and (v\ — v 2 ) — v re i is the relative 
velocity between the cloud and the ambient hot gas. For 
the initial values of setup 2-DD, k m i„ = 0.59 kpc^ 1 which 
corresponds to a wavelength of X m ax = 10.8 kpc. Therefore 
all perturbations with a wavelength shorter than \ ma x are 
leading to KH instability. 

On the other hand, the KH timescale tkh can be esti- 
mated as the time after which the whole cloud mass M gas 
is removed from the cloud by the continuous mass l oss due 
to KH instability. The mass-loss rate is given by iNulsenl 
(|l982l) with Mkh = ^r 2 p am i, x v re i, where r is the cloud 
radius, p am b is the ambient density of the hot halo gas. The 
KH timescale is given by: 



M„ 



trh 



M 



(8) 



KH 



Inserting typical values for our models, the KH 
timescale is 5.8 Gyrs, which is much longer than our simu- 
lation time of 100 Myrs. Starting at 50 kpc with a typical 
velocity of 200 fcms -1 , the cloud reaches the galactic disk 
after 240 Myrs. Therefore, stripping by KH instability does 
not play a major role on the evolution of HVCs. 

2.7.2. Rayleigh Taylor instability 

For an accelerated cloud with a higher density than the 
ambient me dium also the condit ions for a RT instability 
are fulfilled (j Chandrasekhar 196l]). Stabilizing mechanisms 
which can suppress RT instabilities ar e magnetic fields 
(IChandrasekharl["l96ll) . heat conduction (IVieser fc Hensler! 
|2007b|) and self-gravity. iRoediger fc Hensler! (120081 ) investi- 
gated the conditions for RT instabilities in a self-gravitating 
cloud. They found that without magnetic fields, the expres- 
sion for whether a RT instability is present at a radius R 
from the cloud center is equivalent to the expression for 
RPS. Subsequently, in regions where RT instabilities are 
suppressed the cloud is also stable against RPS, while in 
regions with a larger distance to the center, the gravita- 
tional restoring force is smaller which leads to both RPS 
and the occurrence of RT instabilities. 



2.7.3. Velocity dispersion 

As shown in Fig. [5] a significant fraction of the cloud mass is 
located outside the Roche lobe at any time. Therefore, this 



gas is gravitationally unbound to the cloud and very likely 
to be lost within a few dynamical timescales Td yn due to 
the thermal pressure. For a temperature of 1000 K (Setup 
2, 2-D and 2-DD) the mean thermal velocity of the hydro- 



gen atoms within the cloud is an = 5 kr 



and for a 



temperature of 2000 K (Setup 3 and 3-D) <th = 7 kms^ 1 . 
The dynamical timescale, determined by Td yn = D gas /<JH, 
where D gas = 2R gas is therefore between « 50 Myrs (Setup 
2-D) and ~ 105 Myrs (Setup 3). This is of course only valid 
for regions outside the Roche lobe and for a constant tem- 
perature and cloud diameter over time, which is not realis- 
tic. Anyway it provides a first estimate that the simulation 
time of 100 Myrs is of the order of the dynamical timescale 
or twice of it and therefore this process is dominated by 
the ram pressure which is happening on a much shorter 
timescale. 

Unless all other effects mentioned before, this mass-loss 
process is not only valid for the gaseous component but 
also for the DM subhalo that is surroundin g the gas cloud . 
As a typical velocity for the DM minihalo iBurkertl (|1995l) 
finds vq = 17.7(ro/kpc) 2 / 3 km s^ 1 while the virial radius 
as the extensi on of the halo is given by i? v ir = 3.4 x ro 
(|Burkertlll995|) . Subsequently the dynamical timescale for 
the DM component is Td yn ~ 210 Myrs for the ro = 150 pc 
halo and Td yn ~ 250 Myrs for the ro = 250 pc halo. This is 
in both cases more than twice the simulation time, which 
allows us to use a DM halo that is static with respect to the 
cloud center and is therefore accelerated together with the 
gaseous component in the direction of the external mass. 



2.7.4. Tidal disruption 

On their journey to the galactic disk of their host galaxy, 
high velocity clouds would not only be ram pressure 
stripped but also tidally str etched in th e gravi tational field 
of the galaxy. According to iBlitz et all (| 1999T) the already 
disrupted clouds would correspond to the large HVC com- 
plexes with large angular extents and distances in the order 
of a few kpc. A useful estimate f or the cloud's tidal radius is 
provided by the Jacobian limit (jBinnev fc Tremaine!ll988T) 
rj- 



rj = ±D 



-1 1/3 



M(3 + m/M) 



(9) 



where D is the distance between the two mass centers, m 
is the HVC mass and M the galaxy mass. Initially at a dis- 
tance of 50 kpc all cloud radii lie well within the Jacobian 
limit. For Setup 2-DD for example the Jacobian limit is 
4.5 times larger than the cloud radius without a DM sub- 
halo, 10.3 times with a r = 150 pc halo, or 15 times with 
a r = 250 pc halo. After the simulation time of 100 Myrs 
and even with an assumed extreme mass loss of 50% the 
Jacobian limit is still between 2 and 7 times larger than 
the initial cloud radius. In addition, also the virial radius 
of the DM minihalo remains always within the Jacobian 
limit. Therefore we conclude that tidal disruption does not 
play a role on radial orbits and at these distances from the 
galactic center. 
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3. Results 

3.1. Head-tail structure 

For all simulation runs without DM it is possible to re- 
produce the observed and expected head-tail structure in 
the density distribution of the cloud. For a better com- 
parison with observations, the 3D density is integrated in 
z-dircction in order to get a 2D column density distribu- 
tion. Iso-density contours are shown in Fig. [2] as well as 
in Fig. [3] and Fig. 0] While in runs without an additional 
DM potential the iso-density contours are distorted even 
close to the cloud's mass center, they stay very circular in 
models in which the cloud is embedded in a DM subhalo. 
A remarkable difference not only in the head-tail structure 
between simulations with and without DM can be found in 
the results of Setup 2-DD (Fig. [3]) where the cloud without 
DM is totally disrupted while there is hardly any difference 
in the column density distribution between the runs with 
DM halo. We conclude that while the existence of a DM 
halo has a large impact on the head-tail structure of the 
HVCs, it is not very sensitive to the parameters of the DM 
halo as Mdm, and whether the DM density profile has 
a core or a cusp in the center (see Sec. 13.3. 3[) . 

3.2. Line-of-Sight Effects 

Because the HVC motion is by far not radially towards the 
observer, but might happen ballistically, we have chosen dif- 
ferent line-of-sight angles and integrated the radial velocity 
components from a virtual observer. To create a plot in the 
velocity-position space not only the velocity is affected by 
the line-of-sight effect but also the spatial coordinate (see. 
Fig. [5] for further explanations). 

Wc find that different velocity gradients along the head- 
tail structure of the gas cloud can be easily created only by 
changing the inclination of the line-of-sight aspect angle. 
For inclinations of 15deg and 30deg also changes in the 
velocity gradient appear, as observed in different HVCs. 
These changes appear in Setup 3 at distances from the cloud 
center of about 300 pc which is far outside the Roche lobe. 
Therefore the "knees" in phase space according to the dif- 
ferent velocity gradients are not directly connected to the 
Roche lobe, which is at a distance of less than 100 pc (Setup 
3) after 100 Myrs of simulation time. 

For compari son between our models and HVCs obser- 
vatio ns (e.g. bv lBrims et~aTll2000L 120011: iBen Bekhti et al.l 
|2006|) the head-tail structure as a simple horizontal cut 
through the HVCs is not an appropriate measure to decide 
about the existence of a DM subhalo. In the present paper, 
we limit ourselves to only radially infalling HVCs, while 
we suggest that inclined HVC motions - and most prob- 
ably with a significant perpendicular velocity component 
- could lead to the observed knee in the velocity-position 
distribution. 

Since we investigate the location of the HVC in velocity- 
position phase space, we check also the final model for con- 
vergence with runs on higher resolution. The line-of-sight 
dependent contour plots shown in Fig. [5] are analyzed in 
more detail for runs with different refinement levels. The 
mean velocity is calculated in bins of the initial cell size for 
each refinement and shown for an angle of a = 30 deg in 
Fig. [7] for refinement levels 4, 5, 6 and 7. While the gradi- 
ent of the velocity in the tail part is very different between 



refinement runs 4, 5 and 6 the result is nicely converging 
between refinement 6 and 7. 

To compare all simulation runs for their behavior in 
phase space for different line-of-sight directions and inde- 
pendency on their DM content the plots for Setup 2, Setup 
2-D, Setu p 3 and Setup 3-D are shown in Fig. IA.21 Fig. IA.31 
Fig. [O] and Fig.lA~5l 
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Fig. 7. The average line of sight velocity under an inclina- 
tion angle of a — 30 deg in each position bin after 100 Myrs 
of simulation time for four different levels of refinement. 



3.3. Parameter Study 

3.3.1. Effect of the initial wind velocity 

While all models mentioned before start with an initial 
wind velocity of 200 km s _1 i.e. supersonic motion with re- 
spect to the ambient medium of Mach 1.7, in this section, 
we present models of Setup 2-DD without DM halo and 
with a DM halo of rrj = 250 pc at an initial velocity of 
250 /cms -1 (Mach 2.1) in order to investigate the effects 
caused by different relative velocities v re i . 

The gas clouds, however, are not directly experiencing 
the full wind speed due to a quickly forming bow shock 
which leads to a post-shock region with subsonic velocity. 
Fig. [5] demonstrates that only a small difference in the evo- 
lution of the stripping radii exists between the runs with 
the different velocities. 

3.3.2. The dark matter mass 

Setup 2-DD represents a cloud that is in hydrostatic equi- 
librium under its total mass, consisting of the gas {M gas = 
4.4 x 10 5 M ) and a DM minihalo with a scale radius of r$ = 
250 pc and a mass of Mdm = 1-63 x 1O 7 M . Simulations 
with DM minihalos of rn = 150 pc and rrj = 250pc are per- 
formed as well as one model simulation with a pure gas 
cloud only (see Fig. [3J . The maximum refinement level for 
all three runs is set to 7, leading to an effective resolution 
of 5 pc. (For further informations about Setup 2-DD see 
Tab. El) While in the run without DM the cloud is dis- 
rupted and significantly decelerated after not more than 70 
Myrs, both clouds with the Burkert DM halo remain al- 
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most unaffected by ram pressure with only small mass loss. 
They do not show significant differences, even though the 
total DM mass with ro = 250 pc is more than three times 
higher than that of the halo with ro = f 50 pc. 

The time evolution of the model without a DM sub- 
halo is shown by its column density distribution in Fig. IA.1I 
(Appendix). Until 50 Myrs the ram pressure affects only 
the outer HVC shells while the main HVC body remains 
untouched. In this run the central part of the gas cloud is 
highly compressed by means of the strong drag force and 
(starting after 60 Myrs) the peak density is shifted away 
from the simulation rest frame, i.e. is decelerated. 

Due to the very high pressure in the head part of the 
cloud the gas starts to re-expand (perpendicularly to the 
cloud's motion according to the Bernoulli effect) and to 
form a bow-like structure in front of its mass center (be- 
tween 70 and 100 Myrs). This feature shows a compressed 
leading edge but low density in the region behind. Its rela- 
tively large extent makes it very likely to be disrupted by 
both Raylcigh- Taylor instability and RPS in a short time. 

The highest gas mass fraction in our models is applied 
to Setup 3, where the DM mass is only 1.6 times larger 
than the gas mass. Even this light DM halo has a large 
impact on the stability of the cloud. While in the simulation 
run without DM the extent of the HVC perpendicular to 
the wind direction has decreased remarkably within 100 
Myrs, the cloud with the surrounding DM halo conserves 
the initially spherical-symmetric shape much better (see 
Fig. H|) . The mass loss after 100 Myrs is for both runs, with 
and without DM, of the order of a few percent. In Setup 
3 without a DM halo, the ratio between bound mass and 
total mass equals 0.966 which is means that the mass loss 
is twice as large as in the simulation run with a DM halo, 
where the fraction of bound mass is 98.31 %. This means 
that the self-gravity of a gas cloud of 3.1 x 10 6 M Q alone is 
strong enough to keep most of the material bound in the 
first 100 Myrs and cannot be neglected in simulations. 

We determined the bound mass by comparing the sum 
of kinetic and internal energy to the corrected potential 
energy for every cell within the computational domain. 
Because we are interested in the binding energy to the 
cloud or cloud-halo system rather than to the external point 
mass, the total gravitational potential is reduced by the 
potential which is produced by the host galaxy. Every cell 
in which the binding energy dominates contributes to the 
bound mass. 

For Setup 2-DD without a DM halo (Run VI), 95.72% 
of the initial cloud mass is still bound after 100 Myrs. The 
same setup but embedded in an ro = 150pc (Run VII) or 
ro = 250 pc (Run IX) halo reduces the amount of mass loss 
by a factor of 7 - 8 and leads to a bound mass of 99.39% 
(r = 150 pc) and 99.49% (r = 250pc) of the initial gas 
mass after 100 Myrs. 

3.3.3. The dark matter density profile 
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Fig. 8. Semi- analytically derived stripping radius for Setup 
2-DD with DM minihalos of ro = 150pc (full-drawn line, 
Run VI) and ro = 250 pc for an initial relative velocity of 
200 km s^ 1 (dashed line, Run IX) and 250 fcms -1 (dotted 
line, Run XI). 



numerically investigate possible differences in the ram pres- 
sure stripping. 

We re-simulated Run IX, with an NFW halo which 
shows a cuspy density profile unlike the centrally flat 
Burkert halo. We decided to set up an NFW halo with the 
same total DM mass as in Run IX (Af 20 o = 1.63 x 10 7 M ) 
so that we can compare the differences that appear due to 
the shape of the DM halo. The virial radius of this NFW 
halo r2oo = 5.5 kpc and the circular velocity at r2oo is: 
^200 = 3.57 km s _1 . For a full description of the DM density 
profile: 



PNF\v{r) 



wj 1 



(10) 



with S c the overdensity relative to the critical density of the 
universe p c = 3H 2 /(8ttG) (with H = 65 km s' 1 Mpc' 1 ), 
the concentration parameter c in addition to the virial mass 
A/200 is needed. This parameter differs significantly from 
object to object and doesn't show a clear correlation with 
the mass of the object. Fitting rotation curves from galaxies 
from the THINGS survey shows conce ntration parameter 



in large varieties (IChemin et al.l 2011) and vary between 
1 ■ ■ • 30] for dwarf galaxies (|van den Bosch fc Swatersl 



c 

1200% 

Using a concentration parameter of c = 22, which lies 
well within the measured range, leads to a scale radius of 
R s = r 2 oo/c = 0.25 kpc, which is the same as in the com- 
parison halo with a Burkert profile. With the concentration 
parameter c, the overdensity S c is determined by: 



Sr.= 



200 



3 [ln(l + c) - c/(l + c)] 



(11) 



The radial distribution of the DM density as well as the 
circular velocity of both the Burkert and the NFW halo 
are shown in Fig.! 

The results (see Fig. [3]) demonstrate that the final col- 
umn density distribution and therefore the existence of a 
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head-tail structure is neither highly dependent on the DM 
mass (see also Sec. I3.3.2[) nor on the shape of the DM po- 
tential. 
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Fig. 9. Comparison of the DM density profile (bottom) and 
the circular velocity profile (top) for DM halos with a scale 
radius of 250 pc and a mass of 1.63 x 10 7 Af Q . The solid 
line shows the Burkert profile (central core profile), while 
the dashed line represents an NFW profile with a cusp in 
the center. 



4. Conclusions 

Summarizing, the following major issues from the presented 
series of model runs can be derived. While it is possible 
to reproduce the observed change in the velocity gradient, 
this feature is not directly connected to the Roche lobe and 
therefore the position of the observed "knee" in velocity- 
position phase space cannot help to distinguish between 
DM dominated and pure gas clouds. For a variety of initial 
distributions we can create these phase space distributions 
by changing the inclination of the line-of-sight. In a further 
study, clouds with a significant azimuthal motion must be 
tested with respect to their velocity gradient of the stripped 



gas for which we expect a more significantly observable ve- 
locity "knee". 

We simulate DM-embedded clouds with baryonic to to- 
tal mass ratios between 0.03 and 0.39 and found that the 
results depend only marginally on the DM mass and profile 
and that already a DM halo with a mass of 1.6 times the 
baryonic mass stabilizes the gas cloud against ram-pressure 
disruption. 

After 100 Myrs the column-density distributions of the 
clouds with DM halo can be well distinguished from pure 
gas clouds. While the isodensity lines of the core of gas 
clouds without an additional DM gravitational potential 
are losing their initial shape after half the simulation time 
and show a lot of substructures, for all DM-dominated mod- 
els HVCs remain spherical and centrally concentrated, even 
if the DM content is low as in Setup 3. Furthermore an ex- 
tended head-tail structure at galactic heights between 30 
and 50 kpc can only be reproduced for DM-free clouds. 

In our simulations, clouds with DM - with baryonic 
masses in the range from a few times 10 5 to a few times 
10 6 M Q - approach the galactic disk like bullets with high 
velocities, while DM-frce clouds are already decelerated and 
distorted at higher galactic heights. Observed HVCs with 
distances up to 30 kpc, like e.g. Complexes C and H, re- 
veal a very heterogeneous multi-phase internal gas structure 
without any central concentration, which we only obtain in 
simulations without DM. 

Furthermore, the observational lack of HVCs close to 
the Milky Way disk (only LVCs exist) and also that high- 
speed infalling Hi gas is not observed in galaxies with ex- 
tende d Hi envelopes, as e.g. NGC 2403 ([Fraternali et al.l 
I2002D let us conclude that HVCs are decelerated to LVCs 
or disrupted on their infall into galactic halos, contrary to 
the simulated "bullets" of clouds which have an accompa- 
nying DM minihalo. 

Our DM-free models can both explain the proximity of 
the large HVC complexes to the galactic disk even though 
they start their journey from a height of 50 kpc above the 
plane, as well as their morphologies, like e.g. from elongated 
shapes for more massive clouds to head-tail structures and 
the inherent amount of substructures. 

Our conclusions contain the following further serious 
consequences for the cosmological model and for our un- 
derstanding of cosmological structure formation. 

1. Since HVCs are DM-free, they cannot be the cosmolog- 
ical small-scale relics expected from ACDM structure. 

2. Since HVCs have typical distances to the disk of their 
host galaxy of not more than 50 - 60 kpc and do not 
populate the intcrgalactic space frequently enough, it is 
obtruding that, if they are DM-free, they originate from 
stripped-off gas of satellite galaxies or from condensa- 
tions of diffuse intergalactic medium due to thermal in- 
stability rather than from an cosmological background. 
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Fig. A.l. Setup 2-DD. Time evolution of the column density and the equipotential distribution. Contours as in FigO 
From top to bottom, on the left-hand side the panels show snapshots at 0, 20, 30, 40 and 50 Myrs and on the right-hand 
at 60, 70, 80, 90 and 100 Myrs after the simulation onset. 
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Fig. A.3. As in Fig.|X2]but for Setup 2-D. 
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Fig. A.5. As in Fig. [O] but for Setup 3-D. 
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